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Abstract 

We propose an experiment, based on two consecutive Bragg pulses, to mea- 
sure the momentum distribution of quasiparticle excitations in a trapped Bose 
gas at low temperature. With the first pulse one generates a bunch of exci- 
tations carrying momentum q, whose Doppler line is measured by the second 
pulse. We show that this experiment can provide direct access to the ampli- 
tudes u q and v q characterizing the Bogoliubov transformations from particles 
to quasiparticles. We simulate the behavior of the nonuniform gas by numer- 
ically solving the time dependent Gross-Pitaevskii equation. 
PACS numbers: 03.65.-w, 05.30.jp, 32.80.-t, 67.40.Db 
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More than 50 years ago Bogoliubov |IJ developed the microscopic theory of interacting 
Bose gases. A crucial step of the theory is given by the so called Bogoliubov transformations 

bq = Ugdq + V q a\ (1) 

b\ = u q a\ + v q a^ (2) 

which transform particle creation, a, and annihilation, a\ operators into the corresponding 
quasiparticle operators b and W . The real coefficients u q and v q are known as quasiparticle 
amplitudes. The Bogoliubov transformations are the combined effect of gauge symmetry 
breaking and of the interactions which are responsible for the mixing between the particle 
creation and annihilation operators. In virtue of transformations (|l])-(0), the many-body 
Hamiltonian of the interacting Bose gas becomes diagonal in the 6 q 's, representing a system 
of free quasiparticles whose energy is given by the famous Bogoliubov dispersion law: 
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In Eq. (^), c = [gn/m] 1 / 2 is the sound velocity fixed by the density of the gas, n, and 
by the parameter g characterizing the interaction term gY2^( T i ~ r j) °f the many-body 

i<j 

Hamiltonian. The interaction parameter g is determined by the s-wave scattering length a 
through the relation g = Anh 2 a/m. The dispersion law fixes the value of the quasiparticle 
amplitudes u q and v q , which can be written as 

u q ,v q = ±^ ±q2/2m ; (4) 
9 Q 2^t{q) q 2 /2m V ' 

and satisfy the normalization condition u 2 — v 2 = 1. At low momentum transfer (q 2 /2m 

mc 2 ) the Bogoliubov excitations are phonons characterized by the linear dispersion law 

e = qc and the amplitudes u q and v q exhibit the infrared divergence u q ~ —v q ~ (mc/2q) 1 ^ 2 . 

Vice-versa, at high momentum transfer the dispersion law (|3|) approaches the free energy 

q 2 /2m and the Bogoliubov amplitudes take the ideal gas values u q — 1, v q — 0. 

Bogoliubov's theory has been developed also for nonuniform gases. In this case, the 

dispersion law (H) can be defined locally through the density dependence of the sound 



velocity. The theory has been successfully used to interpret the available experimental results 
on the propagation of phonons in trapped Bose-Einstein condensed atomic gases, namely, 
the excitation of the lowest frequency modes 0,0], corresponding to discretized phonon 
oscillations of the system 0,0], the generation of wave packets propagating in the medium 
with the speed of sound || and the excitation of phonons through inelastic photon scattering 
f/|. However, these experiments reveal the propagation of phonons only in coordinate space, 
where the equations of motion take the classical hydrodynamic form, and not in momentum 
space, where Bogoliubov's transformations (]l|)-(H) exhibit their peculiar character. 

In this work we suggest a procedure to measure the Bogoliubov parameters u q and v q in 
a trapped Bose-Einstein condensed gas. Our strategy is based on the following two steps: 

A) First, one generates a bunch of quasiparticles in the sample by means of the technique 
already used in 0. This is based on an inelastic collisional process (two photon Bragg 
scattering) which can be implemented with two detuned lasers transferring momentum q 
and energy Two to the sample. Here q = h(\ti — k 2 ) and uj = {u\ — u^) are fixed by the 
difference of the wave vectors and the corresponding frequencies of the two lasers. In order 
to excite quasiparticles in the phonon regime one should satisfy the condition q < mc. Let 
us call N p h the number of quasiparticles with momentum q generated by this first Bragg 
pulse and let us assume, for simplicity, that the system can be treated as a uniform gas. 
According to the Bogoliubov transformations ([[])- (@), the momentum distribution of the gas 
will be modified as 

n(p) = n (p) + N ph (u 2 q S(p - q) + v 2 q 5{p + q)) , (5) 

where rio(p) is the momentum distribution at equilibrium. Equation (0) reveals the occur- 
rence of two new terms describing particles propagating with directions parallel and an- 
tiparallel to the momentum q of the quasi particles (hereafter called phonons) and weights 
proportional, respectively, to u 2 and v 2 . The total momentum, P = J dp pn(p) carried by 
the system, is equal to qN ph , as a result of the normalization condition u 2 — v 2 = 1. 

B) In the second step of the experiment one measures the momentum distribution (0) 

3 



by sending a second Bragg pulse immediately after the first Bragg pulse. The momentum 
Q, and the energy HQ transferred by the second pulse should be much larger than the ones 
of the first pulse since, in order to be sensitive to the momentum distribution of the sample, 
the scattering should probe the individual motion of particles |8]|| . More precisely, one must 
satisfy the condition hQ ~ Q 2 /2m ^> mc 2 . The measured quantity is the dynamic structure 
factor which, in the large Q regime, takes the form [10| 



Tfi f 

S(Q,Q) = - J dp x dp y n( Px , Py ,p z ) (6) 

where p z = m(hQ — Q 2 /2m)/Q and we have assumed Q to be directed along the z axis. By 
inserting (|5|) into (|]), one finds that the dynamic structure factor exhibits, in addition to 
the original peak located at Q = Q 2 /2mh, two side peaks at 

2mh ma 

By denoting with S + and 5*_ their contributions to the integrated strength J dfl S(Q,Q) = 
N, one finds S + = N ph u 2 and £_ = N ph v 2 or, equivalently, N ph = S + — £_ and v 2 = 
S-/(S + — S-). If the quantity S + + S_ = N ph {u 2 + v 2 ) is much smaller than N, the 
normalization of the central peak remains close to the unperturbed value N. From the 
above discussion one concludes that the measurement of the dynamic structure factor at 
high momentum transfer Q and, in particular, of the two strengths S± would provide direct 
access to the number of phonons generated with the first Bragg pulse, as well as to the value 
of the corresponding quasiparticle amplitudes. 

Expression (^|) for the dynamic structure factor ignores the effects of the final state 
interactions which are responsible for both the line shift of the curve S(Q,Q) and for its 
broadening. These effects can be safely calculated within Bogoliubov's theory and, in the 
large Q domain, are both fixed by the chemical potential of the gas |||9|] . The broadening due 
to mean field effects should not be confused with the Doppler broadening included in Eq. (^j). 
The latter is due to the fact that, even in the equilibrium configuration, the momentum 
distribution of the condensate has a width ~ h/R z originating from its zero point motion 



in the Q direction. In the following, we will consider condensates highly elongated along 
the axial z-axis so that the Doppler broadening, due to the finite size of the system, can 
be ignored. For a safe identification of the two phonon peaks ([7]) and of the corresponding 
strengths S± it is crucial that the separation AQ = ± (q ■ Q) / (mh) between the phonon 
and central peaks be larger than the mean- field effect. This imposes the condition 

qQ/m > fi , (8) 

where we have chosen the two vectors q and Q parallel in order to maximize the separation 
AQ. Equation (H) shows that the momentum q of the phonons generated by the first Bragg 
pulse should not be too small. 

In the second part of the work we explore in detail the microscopic mechanisms of 
generation of phonons produced by the first Bragg pulse, taking into account the fact that 
our system is nonuniform and that the time duration of the pulse is finite. We consider a 
gas of interacting atoms initially confined by a harmonic potential of the form Vh Q (x, y, z) = 
m (lo\(x 2 + y 2 ) + u 2 z z 2 ) /2. The generation of phonons is analyzed through the numerical 
solution of the time dependent Gross-Pitaevskii equation for the order parameter \I/(r, t) || 
in the presence of the additional external potential 

^Braggfo t) = Vf(t) [C0s(qz/h - ut)] (9) 

which reproduces the effects of the inelastic scattering associated with the two photon Bragg 
pulse directed along the axial z direction (see for example Ref. ||11|| ). In Eq. (|9]) the parameter 
V is the strength of the Bragg pulse while the envelope function f(t) was chosen of the 
form f{t) = | [1 + tanh (t/t up )], and f(t) = for t > t B . Here t B is the duration of 
the Bragg pulse, while t up fixes its rise time. By varying the values of q and u of the 
perturbation (Q) different regions of the nonuniform gas are excited with different intensity. 
In fact, according to the Bogoliubov dispersion (3), the condensate is in resonance with 
the periodic perturbation cos(qz/h — ut) for values of the density satisfying the condition 
c 2 = gn/m = [h 2 cu 2 — (q 2 / 2m) 2 ] / q 2 . 
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The ground state, corresponding to the stationary solution of the Gross-Pitaevskii equa- 



tion at large negative times t, was obtained by means of the steepest descent method [12 



For the time dependent solutions we have used a numerical code developed in Ref. flo 
suitable for axially symmetric condensates. The parameter V has been chosen in order to 
generate a number of phonons corresponding to 5 — 10% of the total number of atoms. In 
this way, one produces a visible bunch of excitations whose features can be still described 
using linear response theory. Higher values of V were also considered to explore nonlinear 
effects. The Bragg pulse duration £g was always taken to be significantly less than the os- 
cillation time in the axial direction. This requirement is needed in order to relate the total 
momentum transferred by the photons with the actual momentum carried by the system at 
the end of the first Bragg pulse, thereby ignoring the effects of the external force produced by 
the harmonic potential during the pulse. This condition is well satisfied in the experiment of 
Ref. [0] where the total momentum P z of the condensate was measured after a Bragg pulse. 

An example of the density |\l/(r, t^)| 2 as a function of z and for r± = [x 2 + y 2 ] 1 ^ 2 = 0, 
is shown in Fig. 1. The condensate in this figure has N = 6 x 10 7 sodium atoms confined 
in a trap with ujj_ = 2tt x 150 Hz and uo z = 0.12u±. This corresponds to a Thomas-Fermi 
parameter Na/a± = 10000, where a± = [h/ '(mcjjj] 1 / 2 . We have chosen a duration of the 
Bragg pulse ts = 0.25 x 2n/uj_ (~ 1.7 ms) and intensity V = 1.25huj_. The values of q and 
uj are q = lh/a± and u = 4.13u;_i_. With these parameters we are close to the phonon regime 
(q 2 /2m = 0.02 mc 2 ). 

In Fig. 2 we give the corresponding prediction for the dynamic structure factor S(Q, Q), 
measurable with the second Bragg pulse. This quantity is evaluated in impulse approxima- 
tion and is determined by the longitudinal momentum distribution, as in Eq. (^), 

J dp x dp y n(p x ,p y ,p z ) = f dx'dy'dz'dz e -*P*(*-* , )/ft 

x^*(x',y',z,t B )^(x',y',z',t B ) (10) 

where p z = m(hfl — Q 2 /2m)/Q. Final state interaction effects are ignored in this approx- 
imation, but they do not affect the conclusions of our analysis provided condition @ is 
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satisfied. 

Figure 2 clearly shows the appearance of the two peaks in S(Q, Q) at the frequencies 
predicted by Eq. ([!]). The difference S+ — SL between their strengths gives the number of 
phonons N p h] this turns out to be ~ 6 x 10 6 , i.e., about 10% of the total number of atoms. 
We have verified that the results are independent of the choice of the rise time of the pulse, 
provided t up < 0.05 2tt/uj±. We have also checked that the system responds in a linear 
way, by verifying that the value of P z increases quadratically with V, and, for sufficiently 
long times, the number of phonons generated by the pulse increases linearly with ^ as 
predicted by perturbation theory. Moreover, we point out that condition @, which ensures 
the visibility of the two peaks in S(Q,Q), can be satisfied with reasonable choices of the 
momentum Q of the second Bragg pulse. Taking, for example, the value Q = 21 / um~ 1 || we 
get Qq/m = 36Hlj±, to be compared with the value /i = 25Hlv± of the chemical potential. 
It is finally worth noticing that, since each phonon carries momentum q, their number can 
also be obtained by measuring the total momentum P z after the first Bragg pulse, as done 
in the experiment of Ref. N p h = P z /q. This is useful when q <C mc, since in this case 
S + ~ S- and the difference S + — £_ may be difficult to extract. 

The strengths S± can be used to estimate the value of v^. Our results are given in Fig. 3 
as a function of the first Bragg pulse duration. The three curves have been obtained with 
different choices for the transferred energy and momentum, u and q, but they correspond to 
the same resonant density, i.e., the Bragg pulse excites the system in resonance at the same 
density (~ 0.67 of the central value). The conditions for such resonant behavior have been 
taken from the local density approximation discussed in ||. Due to the different values of 
q, the three curves correspond also to different values of v 2 since this quantity depends on q 
and on the density through the ratio mgn/q 2 , as predicted by Eqs. Our results clearly 

show this effect. In order to make the analysis more quantitative, we also report, for each 
curve, the value predicted by Eq. ([|) with e(q) = Two. In the case of a periodic perturbation 
in a uniform gas the calculated curves of v q should coincide with prediction (g). In our 
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calculations we find that the values of v 2 q exhibit oscillations with frequency 2uo and a slight 
decrease as a function of time. The behavior of v 2 q at short times is the consequence of the 
high frequency components contained in the Fourier transform of the Bragg potential (0), 
whose effects cannot be simply described employing a local density picture. The decrease of 
the signal at larger times is likely the consequence of the diffusion of phonons towards regions 
of lower density as well as of nonlinear effects. Despite these effects Fig. 3 clearly reveals the 
important features predicted by Bogoliubov theory for the quasiparticle amplitude v and, 
in particular, its dependence on the relevant parameters of the system. 

The results of Figs. 1-3 refer to conditions of linear or almost linear regime. It is also 
interesting to explore the response of the condensate to a highly nonlinear perturbation 
generating a number of excitations comparable to the total number of atoms. This can be 
achieved by increasing the strength V of the Bragg pulse. In Fig. 4 we show the dynamic 
structure factor S(Q, Q) calculated after the first Bragg pulse, in conditions of high nonlin- 
earity (V = 25hu)j_). Remarkably, the p z = peak, corresponding to the initial condensate, 
has almost disappeared. In this case, the appearance of additional peaks, associated with 
the second and third harmonics p z = ±2q and p z = ±3g in the longitudinal momentum 
distribution, is clearly visible. 

In conclusion, we have suggested an experimental method to measure Bogoliubov's quasi- 
particle amplitudes in a trapped Bose gas at low temperature. In such an experiment the 
condensate is hit by a sequence of two Bragg pulses. The first (low q momentum trans- 
fer) pulse generates a bunch of phonons which are subsequently mapped in momentum 
space by the second (high Q) pulse. A 3D numerical simulation has allowed us to test our 
predictions and to show that our proposal is compatible with the presently available ex- 
perimental possibilities. This experiment would provide the first direct measurement of the 
Bogoliubov quasiparticle amplitudes, which are of fundamental importance in the theory of 
Bose-Einstein condensation. 
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A.B. would like to warmly thank M. Modugno for providing the code needed for the simula- 
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FIGURES 
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FIG. 2. Dynamic structure factor as a function of p z = m(hQ — Q 2 /2m)/Q after the first Brag 



pulse. 
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FIG. 3. function of time for different choices of q (in units of h/a±) and uj (in units of 

huj±). In order, from top to bottom: q = 1, uj = 4.13; q = 1.5, uj = 6.25 and g = 2, u; = 8.44. The 
values predicted by Eq. (|j) are also reported. 
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FIG. 4. Same as Figure 2. Here V is large in order to drive the system out of the linear regime. 
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